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o 
o 

I Recently the next-to-next-to-leading order QCD corrections and the electro- 

■ weak 0(a) corrections to the Higgs-strahlung processes pp/pp — ► WH/ZH+ 

C3 , X have been calculated. Both types of corrections are of the order of 5-10%. 

In this article the various corrections are briefly discussed and combined into 
state-of-the-art predictions for the cross sections. The theoretical uncertain- 
ties from renormalization/factorization scales and from the parton distribution 
functions are discussed. 

> 
CO 

Q 1. INTRODUCTION 

o , 

At the Tevatron, Higgs-boson production in association with W or Z bosons, pp — > WH/ZH+X, is the 
most promising discovery channel for a SM Higgs particle with a mass below about 135 GeV, where de- 
cays into bb final states are dominant [ QQ. At the pp collider LHC other Higgs-production mechanisms 
play the leading role j3), but nevertheless these Higgs-strahlung processes should be observable. 

Qh! At leading order (LO), the production of a Higgs boson in association with a vector boson, pp — > 

VH + X, (V = W, Z) proceeds through qq annihilation J4|, qq' —> V* — > VH. The next-to-leading 
order (NLO) QCD corrections coincide with those to the Drell-Yan process and increase the cross section 
by about 30% [5|. Beyond NLO, the QCD corrections to VH production differ from those to the Drell- 
Yan process by contributions where the Higgs boson couples to a heavy fermion loop. The impact of 
^ ■ these additional terms is, however, expected to be small in general [6 1. Moreover, for ZH production the 

one-loop-induced process gg — ► ZH contributes at next-to-next-to-leading order (NNLO). The NNLO 
corrections corresponding to the Drell-Yan mechanism as well as the gg — ► ZH contribution have been 
calculated in Ref. [7|. These NNLO corrections further increase the cross section by the order of 5- 
10%. Most important, a successive reduction of the renormalization and factorization scale dependence 
is observed when going from LO to NLO to NNLO. The respective scale uncertainties are about 20% 
(10%), 7% (5%), and 3% (2%) at the Tevatron (LHC). At this level of accuracy, electroweak corrections 
become significant and need to be included to further improve the theoretical prediction. In Ref. [ 8 1 the 
electroweak 0(a) corrections have been calculated; they turn out to be negative and about -5% or -10% 
depending on whether the weak couplings are derived from or a(M|), respectively. In this paper we 
summarize and combine the results of the NNLO corrections of Ref. [7| and of the electroweak 0(a) 
corrections of Ref. [ 8 1. 

The article is organized as follows. In Sects. |^ and we describe the salient features of the QCD 
and electroweak corrections, respectively. Section|^contains explicit numerical results on the corrected 
WH and ZH production cross sections, including a brief discussion of the theoretical uncertainties 
originating from the parton distribution functions (PDFs). Our conclusions are given in Sect.l5~1 



2. QCD CORRECTIONS 

The NNLO corrections, i.e. the contributions at C(a 2 ), to the Drell-Yan process pp/pp — > V* + X 
consist of the following set of radiative corrections: 
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Fig. 1: QCD A"-factors for WH production (i.e. from the sum of W + H and W~ H cross sections) at the LHC (l.h.s.) and 
the Tevatron (r.h.s.). The bands represent the spread of the cross section when the renormalization and factorization scales are 
varied in the range \M V h < Mfl (mf) < 3M V h, the other scale being fixed at ^f(mj?) = M V h- (Taken from Ref. Q.) 



• two-loop corrections to qq — > V*, which have to be multiplied by the Born term, 

• one-loop corrections to the processes qg — > qV* and qq — > gV*, which have to be multiplied by 
the tree-level gq and qq terms, 

• tree-level contributions from qq, qq, qg, gg — > V*+ 2 partons in all possible ways; the sums of 
these diagrams for a given initial and final state have to be squared and added. 

These corrections have been calculated a decade ago in Ref. [9] and have recently been updated [ 10 1. 
They represent a basic building block in the NNLO corrections to VH production. There are, however, 
two other sources of 0{a 2 s ) corrections: 

• irreducible two-loop boxes for qq' — > VH where the Higgs boson couples via heavy-quark loops 
to two gluons that are attached to the q line, 

• the gluon-gluon-initiated mechanism gg — > ZH ITTI at one loop; it is mediated by closed quark 
loops which induce ggZ and ggZH couplings and contributes only to ZH but not to WH pro- 
duction. 

In Ref. |7| the NNLO corrections to VH production have been calculated from the results fTOl on Drell- 
Yan production and completed by the (recalculated) contribution of gg — > ZH. The two-loop contri- 
butions with quark-loop-induced ggZ or ggH couplings are expected to be very small and have been 
neglected. 

The impact of higher-order (HO) QCD corrections is usually quantified by calculating the K- 
factor, which is defined as the ratio between the cross sections for the process at HO (NLO or NNLO), 
with the value of a s and the PDFs evaluated also at HO, and the cross section at LO, with a s and the PDFs 
consistently evaluated also at LO: Kno = (Jho{pP/pP —> VH + X)/aLo(pp/pp — > VH + X). A K- 
factor for the LO cross section, Klo> mav a l so be defined by evaluating the latter at given factorization 
and renormalization scales and normalizing to the LO cross sections evaluated at the central scale, which, 
in our case, is given by hf = = Mvh, where Myn is the invariant mass of the VH system. 

The K-f actors at NLO and NNLO are shown in Fig. Q (solid black lines) for the LHC and the 
Tevatron as a function of the Higgs mass Mjj for the process pp/pp — » WH + X; they are practically 
the same for the process pp/pp — > ZH + X when the contribution of the gg — > ZH component is not 
included. Inclusion of this contribution adds substantially to the uncertainty of the NNLO prediction for 
ZH production. This is because gg — > ZH appears at 0(a 2 ) in LO. 



The scales have been fixed to fj,p = = Myu, and the MRST sets of PDFs for each perturbative 
order (including the NNLO PDFs of Ref. O) 

are used in a consistent manner. 

The NLO iv"-factor is practically constant at the LHC, increasing only from Knlo = 1-27 for 
M H = HO GeV to i^NLo = 1.29 for M H = 300 GeV. The NNLO contributions increase the K-factor 
by a mere 1% for the low Mh value and by 3.5% for the high value. At the Tevatron, the NLO if -factor 
is somewhat higher than at the LHC, enhancing the cross section between i^NLO = 1-35 for Mjj = 110 
GeV and iv~NLO = 1-3 for Mjj = 300 GeV with a monotonic decrease. The NNLO corrections increase 
the ET-factor uniformly by about 10%. Thus, these NNLO corrections are more important at the Tevatron 
than at the LHC. 

The bands around the K-factors represent the cross section uncertainty due to the variation of 
either the renormalization or factorization scale from ^Mvh < ^f{^r) < 3Mvh, with the other 
scale fixed at /xr (/Af) = Myu; the normalization is provided by the production cross section evaluated 
at scales /ip = fiR = Mvh- As can be seen, except from the accidental cancellation of the scale 
dependence of the LO cross section at the LHC, the decrease of the scale variation is strong when going 
from LO to NLO and then to NNLO. For Mjj = 120 GeV, the uncertainty from the scale choice at the 
LHC drops from 10% at LO, to 5% at NLO, and to 2% at NNLO. At the Tevatron and for the same Higgs 
boson mass, the scale uncertainty drops from 20% at LO, to 7% at NLO, and to 3% at NNLO. If this 
variation of the cross section with the two scales is taken as an indication of the uncertainties due to the 
not yet calculated higher-order corrections, one concludes that once the NNLO QCD contributions are 
included in the prediction, the QCD corrections to the cross section for the pp/pp — > VH + X process 
are known at the rather accurate level of 2 to 3% relative to the LO. 

3. ELECTROWEAK CORRECTIONS 

The calculation of the electroweak O(a) corrections, which employs established standard techniques, 
is described in detail in Ref. [8). The virtual one-loop corrections involve a few hundred diagrams, 
including self-energy, vertex, and box corrections. In order to obtain IR-finite corrections, real-photonic 
bremsstrahlung has to be taken into account. In spite of being IR finite, the 0(a) corrections involve 
logarithms of the initial-state quark masses which are due to collinear photon emission. These mass 
singularities are absorbed into the PDFs in exactly the same way as in QCD, viz. by MS factorization. As 
a matter of fact, this requires also the inclusion of the corresponding 0{a) corrections into the DGLAP 
evolution of these distributions and into their fit to experimental data. At present, this full incorporation 
of 0{a) effects in the determination of the quark distributions has not been performed yet. However, 
an approximate inclusion of the 0{a) corrections to the DGLAP evolution shows [ 13 1 that the impact 
of these corrections on the quark distributions in the MS factorization scheme is well below 1%, at least 
in the x range that is relevant for associated VH production at the Tevatron and the LHC. This is also 
supported by a recent analysis of the MRST collaboration [14] who took into account the 0{a) effects 
to the DGLAP equations. 

The size of the 0{a) corrections depends on the employed input-parameter scheme for the cou- 
pling a. This coupling can, for instance, be derived from the fine-structure constant a(0), from the 
effective running QED coupling a(M§ ) at the Z resonance, or from the Fermi constant G M via o>q = 
\/2G ^Myy Sy^ i % ■. The corresponding schemes are known as a(0)-, a(M§ )-, and G^-scheme, respec- 
tively. In contrast to the a(0)-scheme, where the O(a) corrections are sensitive to the non-perturbative 
regime of the hadronic vacuum polarization, in the a(M§)- and G M -schemes these effects are absorbed 
into the coupling constant a. In the G^-scheme large renormalization effects induced by the p-parameter 
are absorbed in addition via ac M - Thus, the G^-scheme is preferable over the two other schemes (at least 
over the a(O)-scheme). 

Figure |2] shows the relative size of the O(a) corrections as a function of the Higgs-boson mass 
for pp — > W + H + X and pp — > ZH + X at the Tevatron. The numerical results have been obtained 
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Fig. 2: Relative electroweak correction 5 as a function of Mh for the total cross section of pp — > W + H + X (l.h.s.) and 
pp — > ZH + X (r.h.s.) at the Tevatron in various input-parameter schemes. (Taken from Ref. 1 8 1.) 



using the CTEQ6L1 fl5ll partem distribution function, but the dependence of the relative electroweak 
correction 5 displayed in Fig.|2]on the PDF is insignificant. Results are presented for the three different 
input-parameter schemes. The corrections in the G M - and a (M§) -schemes are significant and reduce 
the cross section by 5-9% and by 10-15%, respectively. The corrections in the a(0)-scheme differ from 
those in the G^-scheme by 2Ar w 6% and from those in the a{M\ )-scheme by 2Aa(M§) 12%. 
The quantities Ar and Aa(M§) denote, respectively, the radiative corrections to muon decay and the 
correction describing the running of a(Q 2 ) from Q = to Mz (see Ref. [8 1 for details). The fact that the 
relative corrections in the a(0)-scheme are rather small results from accidental cancellations between the 
running of the electromagnetic coupling, which leads to a contribution of about 2Aa(M|) « +12%, 
and other (negative) corrections of non-universal origin. Thus, corrections beyond 0(a) in the q(0)- 
scheme cannot be expected to be suppressed as well. In all schemes, the size of the corrections does not 
depend strongly on the Higgs-boson mass. 

For the LHC the corrections are similar in size to those at the Tevatron and reduce the cross section 
by 5-10% in the G^-scheme and by 12-17% in the a(Mf )-scheme (see Figs. 13 and 14 in Ref. 0). 

In Ref. [ 8 1 the origin of the electroweak corrections was further explored by separating gauge- 
invariant building blocks. It turns out that fermionic contributions (comprising all diagrams with closed 
fermion loops) and remaining bosonic corrections partly compensate each other, but the bosonic cor- 
rections are dominant. The major part of the corrections is of non-universal origin, i.e. the bulk of the 
corrections is not due to coupling modifications, photon radiation, or other universal effects. 

Figure |3] shows the iT-factor after inclusion of both the NNLO QCD and the O(a) electroweak 
corrections for pp/pp — > WH + X and pp/pp — > ZH + X at the Tevatron and the LHC. The larger 
uncertainty band for the ZH production process at the LHC is due to the contribution of gg — ► HZ. 

4. CROSS-SECTION PREDICTIONS 

Figure |4] shows the predictions for the cross sections of WH and ZH production at the LHC and the 
Tevatron, including the NNLO QCD and electroweak 0{a) corrections as discussed in the previous 
sections. At the LHC the process gg — ► ZH adds about 10% to the ZH production cross section, which 
is due to the large gluon flux; at the Tevatron this contribution is negligible. 

Finally, we briefly summarize the discussion [ 8 1 of the uncertainty in the cross-section predictions 
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Fig. 3: if-factors for WH production and ZH production at the LHC (l.h.s.) and the Tevatron (r.h.s.) after inclusion of the 
NNLO QCD and electroweak Q(ot) corrections. Theoretical errors as described in FigureQ 



10 



0.1 



0.01 



. 1 1 1 1 < < 1 1 1 1 1 1 < I 1 




o(pp^VH)[pb] 




: 14 TeV 






; qcd 






QCD+EW 






-• • • . 

gg^ZH 

i . , i , 





1 



0.1 



80 100 120 140 160 180 200 

M H [GeV] 



0.01 



1 1 1 1 1 1 1 1 1 1 1 1 

■ o(pF-^VH)[pb] 
"% Vs"= 1.96 TeV 




\ % WH 






zhV % 

NX NX 
NX NX 
Nx XN 
NX \N 
Xx NX 






NX V 
XX \X 
\s. X.X 
x^ X. 

QCD X, 


S. 

v> 

\s 

\. s 

\\ 




QCD+EW 

r i i 1 t i r 1 t i i 1 i i i 


. . , 1 


\\ 



80 100 120 140 160 180 200 

M H [GeV] 



Fig. 4: Cross-section predictions (in the G M -scheme) for WH and ZH production at the LHC (l.h.s.) and the Tevatron (r.h.s.), 
including NNLO QCD and electroweak O(a) corrections. 



Table 1: Total cross sections (in fb) at the Tevatron (y/s = 1.96 TeV) including NLO QCD and electroweak corrections in 
the G M -scheme for different sets of PDFs. The results include an estimate of the uncertainty due to the parametrization of the 
PDFs as obtained with the CTEQ6 1 15 1 and MRST2001 [ 17 1 eigenvector sets. The renormalization and factorization scales 
have been set to the invariant mass of the Higgs-vector-boson pair, fi — no = Mvh- (Taken from Ref. 1 8 1.) 





pp -> WH + X 


pp^ ZH + X 


M H /GeV 


CTEQ6M [15| 


MRST2001 [17| 


CTEQ6M[15| 


MRST2001 [17| 


100.00 


268.5(1) ± 11 


269.8(1) ± 5.2 


158.9(1) ±6.4 


159.6(1) ± 2.0 


120.00 


143.6(1) ± 6.0 


143.7(1) ± 3.0 


88.20(1) ± 3.6 


88.40(1) ± 1.1 


140.00 


80.92(1) ± 3.5 


80.65(1) ± 1.8 


51.48(1) ±2.1 


51.51(1) ±0.66 


170.00 


36.79(1) ± 1.7 


36.44(1) ±0.91 


24.72(1) ± 1.0 


24.69(1) ± 0.33 


190.00 


22.94(1) ± 1.1 


22.62(1) ± 0.60 


15.73(1) ±0.68 


15.68(1) ±0.21 



Table 2: Same as in TableQ but for the LHC (a/5 = 14 TeV) (Taken from Ref. 0.) 





pp -> WH + X 


pp -> ZH + X 


M H /GeV 


CTEQ6M[15| 


MRST2001 [ 17 1 


CTEQ6M [15| 


MRST2001 [17| 


100.00 


2859(1) ± 96 


2910(1) ± 35 


1539(1) ± 51 


1583(1) ± 19 


120.00 


1633(1) ±55 


1664(1) ± 21 


895(3) ± 30 


9217(3) ± 11 


140.00 


989(3) ± 34 


1010(1) ± 12 


551(2) ± 19 


568.1(2) ±6.7 


170.00 


508(1) ± 18 


519.3(1) ± 6.3 


290(1) ± 10 


299.4(1) ± 3.6 


190.00 


347(1) ± 12 


354.7(2) ± 4.3 


197.8(1) ± 6.9 


204.5(1) ± 2.5 



due to the error in the parametrization of the parton densities (see also [ 16 1). To this end the NLO cross 
section evaluated using the default CTEQ6 [ 15 1 parametrization with the cross section evaluated using 
the MRST2001 [ 17 1 parametrization are compared. The results are collected in TablesHand|2] Both the 
CTEQ and MRST parametrizations include parton-distribution-error packages which provide a quantita- 
tive estimate of the corresponding uncertainties in the cross sections. 1 Using the parton-distribution-error 
packages and comparing the CTEQ and MRST2001 parametrizations, we find that the uncertainty in pre- 
dicting the WH and ZH production processes at the Tevatron and the LHC due to the parametrization 
of the parton densities is less than approximately 5%. 

5. CONCLUSIONS 

After the inclusion of QCD corrections up to NNLO and of the electroweak O(a) corrections, the cross- 
section predictions for WH and ZH production are by now the most precise for Higgs production at 
hadron colliders. The remaining uncertainties should be dominated by renormalization and factorization 
scale dependences and uncertainties in the parton distribution functions, which are of the order of 3% 
and 5%, respectively. These uncertainties may be reduced by forming the ratios of the associated Higgs- 
production cross section with the corresponding Drell-Yan-like W- and Z-boson production channels, i.e. 
by inspecting <?pp/p P ^vH+x I ,,J pp/pp->v+x> rendering their measurements particularly interesting at the 
Tevatron and/or the LHC. 



'in addition, the MRST 1 18 1 parametrization allows to study the uncertainty of the NLO cross section due to the variation of 
a B . For associated WH and ZH hadroproduction, the sensitivity of the theoretical prediction to the variation of a B (a s (M|) = 
0.119 ± 0.02) turns out to be below 2%. 
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